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ABSTRACT 

Aims. We study the connection between spatially resolved star formation and young star clusters across the disc of M 51. 
Methods. We combine star cluster data based on B, V, and /-band Hubble Space Telescope ACS imaging, together with new WFPC2 
(/-band photometry to derive ages, masses, and extinctions of 1580 resolved star clusters using SSP models. This data is combined 
with data on the spatially resolved star formation rates and gas surface densities, as well as Ho- and 20 cm radio-continuum (RC) 
emission, which allows us to study the spatial correlations between star formation and star clusters. Two-point autocorrelation func- 
tions are used to study the clustering of star clusters as a function of spatial scale and age. 

Results. We find that the clustering of star clusters among themselves decreases both with spatial scale and age, consistent with 
hierarchical star formation. The slope of the autocorrelation functions are consistent with projected fractal dimensions in the range of 
1.2-1.6, which is similar to other galaxies, therefore suggesting that the fractal dimension of hierarchical star formation is universal. 
Both star and cluster formation peak at a galactocentric radius of ~2.5 and ~5 kpc, which we tentatively attribute to the presence of 
the 4:1 resonance and the co-rotation radius. The positions of the youngest (< 10 Myr) star clusters show the strongest correlation 
with the spiral arms, Ha, and the RC emission, and these correlations decrease with age. The azimuthal distribution of clusters in 
terms of kinematic age away from the spiral arms indicates that the majority of the clusters formed ~ 5-20 Myr before their parental 
gas cloud reached the centre of the spiral arm. 

Key words, galaxies: individual: M51 - galaxies: star clusters 



1. Introduction 

The increase in the amount of multi-wavelength data during 
the past few decades has been of tremendous importance for 
the field of star formation. On the largest scales, these multi- 
wavelength studies focussed primarily on the tight correlation 
between global radio-conti nuum (RC) and far-infrared (FIR) 
emission of galaxies (e . g. | Condon|[l992l : iNiklas & Beck|[l997t 
1 lNiklas|[l997riYun et alj[200ll) and on the global Schmidt-law, 
which dictates that there is a powerlaw relation between the 
area-normalised star formation rate (SFR) and total gas sur - 
face density with an index of 1.4 + 0.15 dKennicuttl Tl998b). 
These correlations have been proven useful in the derivation of 
glob al star formatio n rates from radio-continuum observations 
(e.g. ICondonlll992l). as w ell as in models of galaxy evolution 
(Mih os & Hernquisllll994l) . 

More recently, these multi-wavelength correlations are ex- 
ploited to study star formation locally across single galax- 
ies. This is done by either utilizing the RC-FIR correlatio n 
dMurphv et al.ll2006l: ISchuster et aLll2007t iMurphv et al]|2008l) . 
or by combining Spitzer Space Telescope FIR observations with 
Ha and Paa data, which allows for a more tho rough extinc- 
tion correction on the derive d star formation rates dCalzetti et alj 
120051: iKennicuttetal] 120071) . Combined with new CO and HI 
data, recent studies show that the Schmi dt-law is also valid 
locally across the s piral galaxy M 51 (Schuster et al. 2007; 
iKennicutt et~aT]|2007l) . 

There is ample evidence that star f ormation leads to forma- 
tion of star clusters (Lada & Lada 2003, and references therein). 
Therefore, star formation across spiral galaxies has been stud- 



ied extensively through their star cluster populations , eithe r 
by using ground-based data dLarsen & Richtlerl 1 19991 |2000|) . 
or by exploiting the highe r resolution Hubble Space Telescope 
(HST) observations (e.g. Holtzman et al. 19921: Meurer etail 
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119951: IWhitmore et aljH999t Ide Griis et aTll2003ULarsenll2004 
Most of these studies, however, focus on the star cluster prop- 
erties themselves, on the connection between the galaxy's star 
cluster populati on and SFR on gl obal scales, or on their spa- 
tial correlation dZhang et al.l l2001). There is still a lack of evi- 
dence, however, on the actual fraction of cu rrent star formation 
taking place in young star clusters (although lde Griis et al. i 2003 
find a lower limit of « 35% for the Mice galaxies and Fal l et alj 
2005 find a lower limit of 20% in the Antennae). Although it 
is probably true that most star fo rmation takes place in some 
sort of embedded, clustered form dLada & L ada 2003), it is not 
clear yet which fraction of these clusters are actually observ- 
able as star clusters after they emerge from their natal gas cloud 
dBastian & Goodwinll2006l:lGoodwin & Bastianll2006l) and if, or 
how, this depends on environment. The answers to such ques- 
tions could give us important clues on the process of star forma- 
tion and the early evolution of young star clusters. 

Before such a quantitative comparison can be made, how- 
ever, we first need to determine the amount of spatial corre- 
lation between star and cluster formation, since any quantita- 
tive comparison is only meaningful for regions that are spatially 
correlated. Besides, determining the spatial correlation between 
clusters themselves allows us to study the amount and evolu- 
tion of clustering of star clusters, which provide important clues 
as to which extent the formation of clusters is hi erarchical (e.g. 
lEfremov & Elmegreenlll998UBastian et~al"1l2007l) . 

For M 5 1 data with a wide wavelength range and good spatial 
resolution are available. This makes it possible to study, first of 
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all, how well the positions of star clusters are tracing the regions 
of star formation across this spiral galaxy, and second of all, 
which fraction of the star formation is taking place in the form of 
visually observable star clusters. In this study we combine data 
on th e resolved Schmidt-law in M 5 1 (SFR and gas surface den- 
sities, |KennicuttetaniQ07l) as well as archival radio-continuum, 
2MASS and HST/WFPC2 data with slig htly resolved young star 
clusters from lScheepmaker et al.l 1 120071) . These datasets are cov- 
ering a large fraction of the disc, allowing us to compare star and 
cluster formation under a range of conditions, such a low/high 
gas and stellar surface densities, within a single galaxy consis- 
tently. 

In this work we will introduce the new datasets used, and 
present new WFPC2 Uf336w photometry. This photometry al- 
lows us to add ages, masses and ex tinctions to the star clus- 
ter data of Scheepmake r et al] d2007l) . In the remainder of this 
work we will then focus on the spatial correlations between star 
and cluster formation, for differently aged cluster populations. 
A quantitative comparison between star and cluster formation, 
as well as the cluster initial mass function (CIMF), w ill be the 
subject of a follow-up paper (IScheepmaker e t al. 2009). 

In Sect. |2] we describe the star cluster data we used, as well 
as the HST/WFPC2 observations. In Sect. [3] we estimate ages, 
masses and extinctions of the star clusters in order to select our 
cluster samples. We describe the data on SFRs and gas densities 
in Sect. [4] together with our method of tracing the spiral arms. 
Finally, in Sect.[5]we study the spatial distribution of star clusters 
and star formation in M 5 1 and we summarize our main results 
in Sect.|6] 




Right Ascention (J2000) 



2. Observations and photometry 

For our star cluster data, we make use of the dataset of 
Scheepmaker et al. (2007, from here on referred to as "S07"), 
complemented with HST/WFPC2 U-band (Uf336w) observa- 
tions. 



2.1. HST/ACS observations 

S07 used the HST/ACS mosaic image of M 51 in B F435W (~B), 
Vfsssw (~V) and Ifs\4w (~-D to select star clusters across the 
complete disc of M 51, based on size measurements. For a 
comple te description we refe r to S07 for the star cluster data 
and to Mutch leret al.l d2005l) and the M 51 mosaic website 
(http://archive.stsci.edu/prepds/m51/) for the image 
processing. Below, the source selection, photometry and size 

measurements of the clusters are summarized. 

The SExtractor package dBertin & Arnoutsl 1 19961) was 
used to select sources above a local threshold, cross- 
correlating between the three filters. Aperture photometry in 
the VEGAmag system was performed on the sources using the 
IRAFlJ'DAOPHOT package, using an aperture radius of 5 pix- 
els and a background annulus with an inner radius of 10 pix- 
els and a width of 3 pixels. A fixed aperture correction for a 
3 pc source was applied, which was respectively -0.17 mag for 
Bf435w and Vf555w and -0.19 mag for If&uw, and the photome- 
try was corrected for Galactic foreground extinction in the direc - 
tion of M 5 1 according to Appendix B of ISchlegel et all (Il998l) . 
The star clusters, being slightly extended, were distinguished 



Fig. 1. Footprints of the WFPC2 Uf336w fields used in this study, 
overlaid on the ACS Bp435w mosaic image. The numbers refer to 
the different datasets in Table Q] 



from stars based on siz e measurements using the Ishape rou- 
tine dLarsenll 19991 120041) . assuming EFF 15 models dElson et all 
1 1987b for the surface brightnesss profiles of the clusters. 



1 The Image Redcution and Analysis Facility (IRAF) is distributed by 
the National Optical Astronomy Observatories, which are operated by 
the Association of Universities for Research in Astronomy, Inc., under 
cooperative agreement with the National Science Foundation. 



2.2. HST/WFPC2 observations 

For a reliable estimation of star cluster para meters using Simple 
Stellar Population models (see Sect. |3~1| >, And ers et al.l ((2004) 
have shown that at least four passbands are necessary, with the U 
and B bands having the highest significance for young clusters. 
We therefore complemented the HST/ACS data with archival 
HST/WFPC2 observations in the Uf336w passband (~ U). 

A total of seven pointings were taken from the ESO/ST- 
ECF science archive. The footprints of the pointings, overlaid 
on the ACS fif435w image, are shown in Fig. Q] An overview 
of the different datasets is given in Table Q] We note that pro- 
posal ID 10501 includes another Uf336w pointing covering the 
companion galaxy of M 51, NGC 5195. Since we are interested 
in the star clusters in the disc of M 51, this pointing (dataset 
U9GA0101B) was not used. We also did not use the WF4 chip 
of dataset U9GA0601B (field 7 in FigQ}, since it contained too 
few sources. This made the manual selection of enough refer- 
ence objects impossible, which were necessary to transform the 
cluster coordinates from the ACS frame to the WFPC2 frames 
(see below). 

We used the individual chips (WF and PC) of the standard 
pipeline-reduced associations, i.e. no further data reduction was 
performed. On every chip between 6-1 1 sources were matched 
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Field" 


ID 


PI 


Dataset (association) 


Exposure time 


Observation date 


1 


5652 


R. R Kirshner 


U2EW0205B 


3 X400 = 1200 s 


12 May 1994 


2 


7375 


N. Scoville 


U50Q0105B 


2 x 600 = 1200 s 


21 Jul 1999 


3 






U9GA0201B 


2 X 1300 = 2600 s 


24 Jan 2006 


4 






U9GA5301B 


2 x 1200 = 2400 s 


28 May 2006 


5 


10501 


R. Chandar 


U9GA0401B 


2 x 1300 = 2600 s 


13 Nov 2005 


6 






U9GA0501B 


2 x 1300 = 2600 s 


13 Nov 2005 


7 






U9GA0601B 


2x 1300 = 2600 s 


13 Nov 2005 



" The numbers refer to the footprints in Fig. |JJ 



with sources on the ACS Bf435w image. Using these sources the 
coordinate transformations between the WFPC2 chips and the 
ACS frame were calculated using the IRAF task geomap. The 
task geoxytran was used to transform the coordinates of the 
7698 resolved sources of S07 on the ACS Bfmsw image to the 
WFPC2 chips. In most cases the root-mean-square accuracy in 
X and Y was better than one pixel, and in all cases better than 2 
pixels. 

Aperture photometry was performed at the transformed co- 
ordinates using DAOPHOT, allowing for a maximum centring 
shift of three pixels. For the PC chips a (5, 10, 3) configuration 
was used for the radius of respectively the aperture, annulus 
and width of the annulus in pixels. For the WF chips the con- 
figuration was (3,7,3). Aperture corrections to an aperture of 
0.5"were measured using analytic cluster profiles with an effec- 
tive radius of 3 pc, which were subsequently convolved with the 
PSF These corrections for the WF and PC chips were -0.14 and 
-0.19 mag for the WF and PC chips, respectively. Zeropoints 
(calibrated for the 0. 5"aperture) and CTE loss corrections were 
calculated following lDolphin l ( 2000 ) We corrected for Galactic 
foreground ex tinction in the dire ction of M 51 according to 
Appendix B of lSchlegel etai] dl998l) . 

From the initial source list of 7698 resolved sources on the 
ACS Bf435w image, we retrieved positive photometry on the 
WFPC2 Uf3j,(,w images of 5502 sources. The remaining sources 
had either an aperture covering the edge of an image, a failure 
in daophot' s centring algorithm or a negative flux, caused by a 
combination of a faint source and noise in the background an- 
nulus. As can be seen in Fig. Q] there is some overlap between 
different WFPC2 chips. If a source had photometry in more than 
one Uf33(,w chip, the photometry with the smallest error was se- 
lected. 

In Fig. [2] we show the resulting photometry and error in the 
photometry for the Uf336w passband. As expected, the photome- 
try of fields 1 and 2 has larger errors compared to the photometry 
of fields 3-7, due to the deeper exposures of the latter fields. The 
horizontal line at ct(Uf336w) = 0.2 mag indicates the limit we 
applied to select only clusters with accurate photometry. The bin 
centred on Uf336w = 22 mag indicates the magnitude at which 
approximately 10% of the clusters will be rejected from the sam- 
ple by this ct(Uf33(,w) < 0.2 mag criterion. These limits will be 
used in Sect. !3.3l to select cluster samples. 

3. Selection of the cluster samples 

In this study we are primarily interested in the spatial distribu- 
tion of young star clusters in the disc of M 51. We therefore 

2 The most recent equations for the CTE losses were taken from 
http : //purcell . as . arizona . edu/wfpc2_calib/. 
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Fig. 2. The error in the £/f336w magnitude versus the Uf336w 
magnitude for all 5502 sources for which we have photome- 
try. The different symbols/colours refer to the different datasets 
listed in Table Q] The horizontal line indicates ctUf336w = 
0.2 mag and the bin at Uf336w — 22 mag indicates the magnitude 
range in which ~90% of the clusters has ctUf336w < 0.2 mag. 

used Simple Stellar Population (SSP) models to estimate the age, 
mass and extinction of all 5502 clusters for which we have UBVI 
photometry (Sect. [XT). The estimated ages were then used to se- 
lect cluster samples with different ages (Sect. [33). 

3.1. SSP model fitting 

For our age determinations we used the galev SSP m odels 
(ISchulz et all l2002t lAnders & Fritze-v. Alyenslevenl I2003I) to- 
gether with the analysed tool ( I Anders et al.ll2.004l) . In summary, 
analysed compares the broad-band spectral energy distribution 
(SED, in our case UBVI) of every cluster to a large grid of 
modelled SEDs for different ages, extinctions and metallici- 
ties. The model with the lowest x 1 is selected, and the (initial) 
mass of the cluster is derived by scaling the entire modelled 
SED to the observed S EP. The models and the tool have been 
extensively tested (e.g. lAnders et al.ll2004t Ide Griis et aL1l2005b 
Ide Griis & And ers 2006) and the uncertainty in the derived ages 
is typically < 0.4 dex. 

The galev models were preferred above other SSP mod- 
els because of their freedom in stellar evolutionary models and 
metallicity and because their output and treatment of extinc- 
tion is tailored to our specific HST/ACS filters. Therefore, no 
transformation between different photometric systems was nec- 
essary. We used the models based on the Geneva isochrones 
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dCharbonnel et alj|1993t ISchaerer et alj|1993l) . which are better 
suited for the age range we are interested in (i.e. age < 100 Myr) 
than e.g. Padova isochrones dBertelli et al.l Il994t iGirardi et al.l 
2000), due to a higher time resolution in this age range. However, 
in Sect. l3.2l we will discuss the impact of the different isochrones 
on the derived cluster ages. For the models we further adopted a 
lKroupal d2001) stellar initial mass function (IMF) down to 0. 1M 
and a metallicity between 0.4Z o and 2.5Z . This range in metal- 
licity covers the observed met allicity of H II regions in the disc 
of M 51 (solar to twice solar. iDiaz et~aLlll99 It Hill et al.ll 1997b 
and the slightly lower metallicity expected for older clusters. The 
extinction w as varied between < E(B-V) < 1 in steps of 
0.05 and the ICardelli et alj d!989l) extinction law was applied. 
Absolute magnitudes of all the clusters were calculated assum- 
ing m - M - 29.62 ( i.e. a distance to M 51 of 8.4 + 0.6 Mpc, 
Feldmeie r et al.|[l997h . In Fig.[3]we compare the photometry of 
our clusters to the SSP models in two colour-colour diagrams. 



3.2. Systematic uncertainties 

Due to either the discreteness of the model isochrones or the 
rapid colour evolution of the models at certain ages in com- 
bination with the photometric uncertainties of the data, some 
fitting artefacts, or "chimney's" will gen erally appear in the 
log(mass) vs. log(a ge) plane of the clusters (Bas tian et al.l20 05a; 
Giel es et ai1l2005l) . We find that these artefacts depend on the 
SSP model parameters such as the adopted isochrones and the 
range in metallicity. In Fig. [4] we compare the ages of clusters 
fitted with our adopted parameters (Geneva isochrones, Z free 
between O.4Z - 2.5Z Q ) to the ages fitted with metallicity fixed 
to solar or using Padova isochrones. 

First of all, Fig. [4] (left) shows that restricting the metallic- 
ity to be solar, leads to a more discrete distribution of ages and 
the exchange of clusters between the different chimney's: many 
clusters from the log(age/yr) = 7.2 chimney (on the x-axis) are 
now fitted with log(age/yr) between 7.6-8 or around 6.7 (on the 
y-axis), i.e. they show a vertical scatter. Similarly, many clus- 
ters from the range log(age/yr) = 6.7-8.0 (x-axis) end up in the 
log(age/yr) = 7.2 chimney (y-axis), representing the horizontal 
scatter. 

Secondly, Fig.|4](middle & right) shows that with the Padova 
isochrones we generally fit older cluster ages for log(age) $7.5 
(x-axis). These panels also show the limited time resolution 
of the galev models based on Padova isochrones, in the first 
~ 16 Myr. This is caused by the fact that in the galev models the 
resolution can only be as high as the lowest possible age, which 
is 4 Myr when using the Padova isochrones. The preference of 
Geneva based models for younger ages is likely a combination 
of a stronger contribution from red supergiants, leading to red- 
der colours, and the higher age resolution, which lead to less 
"smoothing" of the colour fluctuations compared to the Padova 
based models. 

We conclude from Fig. [4] that cluster age determinations 
based on UB VI broad-band photometry can be very uncertain in 
the first ~ 100 Myr due to the uncertainties in the adopted model 
parameters. These uncertainties will translate into uncertainties 
in the different age samples, which we will select in Sect. 13.31 
Although we will use the ages based on the models using Geneva 
isochrones and O.4Z < Z < 2.5Z Q as our "standard ages" in the 
remainder of this work, we will check, when necessary, how ro- 
bust our results are against using the other models of Fig. |4] 



3.3. Selection of the different age samples 

After fitting the SSP models to our sample of 5502 clusters with 
UBVI photometry, we applied the following selection criteria to 
reject clusters with unreliable fits and to minimize incomplete- 
ness effects: 

1. Photometry in Bf^sw, Vfsssw and Ifsuw brighter than the 
90% magnitude completeness limits determined by S07 for 
a source with an effective (half-light) radius of 3 pc in a high 
background region (respectively 24.2, 23.8 and 22.7 mag). 

2. Photometric accuracy in Bp43sw, Vfsssw, Ifsuw an d Uf336w 
better than 0.2 mag. 

These two criteria reduced our cluster sample from 5502 to 
2089 clusters. Since the Uf336w imaging was not part of the 
source detection, we do not have a well defined completeness 
limit in the Uf336w passband. We can, however, make a rough 
estimate of the 90% completeness limit in U F336W, since we per- 
formed the U F33(,w aperture photometry at the coordinates of all 
the resolved clusters on the fi/?435w image. This photometry was 
followed by the (tUf336w < 0.2 mag criterion. Using bins of 
0.2 mag wide, we measured that the magnitude at which this 
(tUf336w criterion recovers ~9Q% of the clusters from the to- 
tal sample is Uf336w - 22.0 (see also Fig. |2j. This magnitude 
can be considered as a rough estimate of the U F336W 90% com- 
pleteness limit, since in a real completeness test the recovered 
fraction is measured as a function of the input magnitude and 
not of the measured (possibly extincted) magnitude. Also, in a 
proper completeness test the recovered clusters are compared to 
an input sample which is (by definition) 100% complete. In our 
case we can only measure the effect of the <tUf33(,w criterion 
on the detected sample, which is probably already affected by 
some incompleteness due to extinction. Keeping this in mind, 
we used the estimated U F33W limit as an extra criterion to select 
our cluster samples: 

3- Uf33(,w < 22.0 mag, 

which further reduced our sample to 1580 clusters. 

To estimate the impact of an extinction induced bias, we can 
do the following: if we would assume that Uf336w = 22.0 mag 
is the detection limit for extinction free sources, we could cor- 
rect the measured Uf336w photometry of every source accord- 
ing to their best-fitting extinction, before applying the magni- 
tude cut. This would include another 271 clusters in our sample. 
Still, such a correction would not account for the most extincted 
sources, which are not detected in the first place. This shows that 
extinction effects can introduce biases of the order of ~ 15% of 
the sample size. Unfortunately we have to accept that we sim- 
ply can not fully correct for such extinction effects by the proper 
choice of magnitude cuts. 

The extinction distribution is sho wn in Fig. [5] an d is in ex- 
cellent agreement with the results of lBik et all (12003). Since the 
distribution is peaked towards low extinctions, most of the ex- 
tra clusters, if we could correct for extinction, would be located 
close to the Uf33w detection limit. We therefore note that, al- 
though one can use 90% completeness limits, since extinction 
(and also crowding, see e.g. S07) can not properly be taken into 
account, one will always end up with a sample being less than 
90% complete. This also holds for the remainder of this work, 
where we nevertheless will use terms like "90% complete" for 
the sake of simplicity. 

The age-mass diagram for the clusters fullfilling criteria 1-3, 
is shown in Fig. [6] In this diagram two trends are immediately 
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Fig. 3. U F336W - B F4 35w versus V F555W - Ipsuw diagram (left) and B F435W - V F555W versus V F555W - I F guw diagram (right) of the 
1580 clusters selected in Sect. l3.3l (not corrected for extinction), compared to galev SSP models using Geneva isochrones, a Kroupa 
IMF and three different metallicities, indicated in the figures. The stars indicated log(age/yr), starting from 6.6 (left) to 9.8 (top 
right) in steps of 0.4 dex. 
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Fig. 4. Comparison of cluster ages determined using different galev SSP model parameters, for the 1580 clusters selected in 
Sect. 13.31 The line shows the one-to-one correlation. Left: Geneva isochrones for solar metallicity versus Geneva isochrones for 
Z < Z < 2.5Z . Middle: Padova isochrones versus Geneva isochrones, both for Z < Z < 2.5Z Q . Right: Padova isochrones for 
solar metallicity versus Geneva isochrones forZ < Z < 2.5Z . 



visible. First, we see that the mass of the most massive clus- 
ter generally increases with age. This is most likely a statisti- 
cal size- of-sample effect caused by the logarithmic age interval 
(see e.g. lHunter et alj2003l:lGieles & Bastianl2008l) . Second, we 
see that the mass of the least massive clusters increases steadily 
with age. This is a selection effect caused by the detection limit 
(i.e. the most limiting 90% magnitude completeness limit, see 
criterium 1 above). For the BVI passbands we also plotted the 
mass evolution with age, corresponding to the 90% magnitude 
completeness limit. For the U F336W passband we plotted the cor- 
responding mass evolution for the applied U F33 (,w = 22.0 limit. 
We see that for log(age/yr)> 6.7 U F33 (,w is the most limiting fac- 
tor that sets the lower mass limit. These mass evolution curves 
were calculated using galev models with solar metallicity. The 
slope of these detection limits is caused by the evolutionary fad- 
ing of the clusters: a less massive cluster reaches a constant mag- 
nitude limit much faster compared to a more massive cluster. 
These curves show that after ~ 20 Myr the limiting mass in 
U F33 w evolves faster with age compared to the other passbands, 
because the integrated flux in U F33 (,w is always dominated by the 
hottest (and therefore fastest evolving) stars. We see some clus- 



ters with best-fitted masses below the predicted limiting U F33 (,w 
mass. This is possible because of two effects. First of all, the 
mass of a cluster is determined by scaling of the complete mod- 
elled SED (i.e. UBVI) and not by scaling a single passband only. 
Secondly, the curve for the limiting mass is for solar metallicity, 
while the best-fitting model can have a different metallicity. 

We selected three magnitude -limited cluster samples with 
different ages for the remainder of this study: 

• Sample 1: 

6.4 < log(age/yr) < 7.0 
Number of clusters: 895 

• Sample 2: 

7.0 < log(age/yr) < 7.5 
Number of clusters: 556 

• Sample 3: 

7.5 < log(age/yr) < 8.6 
Number of clusters: 127 

These three cluster samples are indicated by the numbered re- 
gions in Fig. [6] The total number of clusters in these three sam- 
ples is 1578. 
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Fig. 6. Log(mass) versus log(age) diagram of the 1580 clus- 
ters with photometric accuracies < 0.2 mag and extinction cor- 
rected magnitudes brighter than the 90% completeness limits in 
Bf435w, Vf555w and Ifuaw and Ufh6w < 22.0, as indicated in 
the lower-right corner. The corresponding mass limits as a func- 
tion of age are indicated by the coloured lines. Note that the 
masses are initial masses (i.e. corrected for stellar evolution but 
not for secular evolution). The three magnitude-limited age sam- 
ples are indicated by the numbered regions. 



Due to the systematic effects described in Sect. 13.21 our clus- 
ter samples will change significantly if we adopt the Padova stel- 
lar isochrones instead of the Geneva isochrones. This results in 
~ 230 clusters less in sample 1, and ~ 200 clusters more in sam- 
ple 3. 

As can be seen from Fig. [6] our cluster samples have different 
lower-mass limits. Our youngest cluster sample contains clus- 
ters down to estimated masses of log(M/M )w 2.6, but is 90% 
complete down to log(M/M ) = 3.1, while the lowest possible 
mass in our oldest cluster sample is log(M '/ M G )^ 3.7. An elec- 
tronic table with the photometry, ages, masses and extinctions of 
the clusters will be publis hed together with a forthcoming paper 
dScheepmaker et al.ll2009l) . 



4. Star formation rates, gas densities and spiral 
arms 

4. 1 . Spiral arms 

In order to study the relation between star and cluster formation 
and spiral structure, we first defined the two spiral arms in M 5 1 
using data from the Two Micron All Sky Survey (2MASS), which 
we downloaded from the NASA/IPAC Extragalactic Database 
(NED). The background subtracted //-band (1.65/mi) image 
was used, which is part of the 2MASS Large Galaxy Atlas 
(Jarre tt et al. 1 120031) . We used near-infrared data to trace the spi- 
ral arms, since the near-infrared is most sensitive to older stellar 
populations. This means that we are tracing the dominant mass 
component in the disc of M 51, independent of the dominant 
flux in the near-UV and optical of young star clusters and mas- 
sive star formation. 

We enhanced the spiral structure by first masking the com- 
panion galaxy, the inner 750 pc region, bright point sources and 
negative pixels. We then subtracted radial averages and con- 
volved the resulting image with a Gaussian kernel with a FWHM 
of 25" (* 1 kpc). In the resulting smoothed image we identified 
23 peaks tracing the two spiral arms using IRAF/DAOPHOT, 
and transformed their positions into polar coordinates. Finally, 
the regions in between the peaks were interpolated using cu- 
bic splines of the form ln(R(ff)). The resulting interpolated spiral 
arms are shown in Fig.[7j overplotted on the original background 
subtracted H-band image (left) and the enhanced image (right). 
We will refer to the arm extending towards the south-west (i.e. 
bottom-right) as "Spiral arm 1", and to the arm extending to- 
wards the nort-east (top-left) as "Spiral arm 2", as indicated in 
the figure. 

4.2. Star formation rates and gas surface densities 

In a recent paper lKennicutt et"ai] (|2007, from here on referred to 
as "K07") studied the locally resolved SFRs and gas densities in 
M 51, based on Ha, Paa, CO, HI and Spitzer Space Telescope 
24pm imaging. For a series of 257 circular apertures with a di- 
ameter of 13" (corresponding to ~530 pc) and distributed across 
the disc of M 51, these authors present extinction corrected SFR 
surface densities, as well as hydrogen gas and total gas surface 
densities. The SFRs are derived from the Ha line emission (after 
extinction correct ion using Paa) using the transformation from 
iKennicutl dl998al) . and a re therefore tracing the current star for- 
mation up to -10 Myr dCalzetti et alj|2005l) . We used Table 1 
of K07 to retrieve their list of aperture positions and the corre- 
sponding SFR surface densities and total gas surface densities. 

The aperture positions were transformed to the frame of the 
ACS Vf555w image using the IRAF task rd2xy. In Fig. [8] (left) 
we show the resulting positions of the 13" apertures. The mean 
SFR surface density in every aperture is indicated in greyscale. 
In the regions where two or more apertures are overlapping the 
average value was used. 



5. Spatial distribution of star clusters and star 
formation 

In Fig. [8] (right) we show the 2-dimensional spatial distribution 
of the three cluster samples in M 5 1 . Also shown are the interpo- 
lated spiral arms from Fig. [7] and the regions (white) which are 
covered by the WFPC2 pointings, i.e. the regions for which we 
have derived ages of the clusters. 
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Fig. 7. Left: Location of the two spiral arms, overlaid on the 2MASS Large Galaxy Atlas H-band image of M 51. The few black 
points are foreground stars. Right: Same spiral arms, overlaid on the enhanced (masked, radial average subtracted and smoothed) 
image. 



We can see from Fig. [8] that most current star formation is 
taking place in the centre of M 5 1 and at a galactocentric radius 
of ~ 5 kpc, and that these are also the locations where most of 
the young star clusters reside. We study these correlations more 
closely in Sect. 15. II with the help of radial distributions. 

Fig. [8] shows clearly that the majority of the clusters is 
closely tracing the spiral arms. This is expected, since the ma- 
jority of the clusters are believed to be formed by the passing 
spiral density wave. Older clusters are expected to trace the spi- 
ral arms less closely, and we indeed see that the clusters from 
sample 3 (log(age/yr) > 7.5) are more uniformly distri buted. 
This is consistent with the result of Hwan g"& Led ([2008), who 
found that, using only the ACS BVI data, the bluest clusters in 
M 5 1 are closely associated with the spiral structure and redder 
clusters are more uniformly distributed. We look at this in more 
detail in Sect. l5.2l using azimuthal distance distributions. 



East (left) of the centre of M 5 1 the clusters are more con- 
centrated on the convex (outward) side of spiral arm 1 compared 
to the concave (inward) side, up to ~5 kpc. For spiral arm 2 in 
the south-west this effect is less obvious, although we note that 
for spiral arm 2 we are much more hampered by selection effects 
due to the incomplete WFPC2 coverage. 

Fig. [8] also shows that the distribution of clusters it- 
self is "clustered" in so-called cluste r complexes (see 
e.g. lEfremovl [1995 : Bastian et alj l2005bl) . These complexes 
are the largest groups in the hierarchy of star formation, 
whic h originates from the fractal distribution of the interstellar 
gas dElmegreen & Efremovll 19961 : lEfremov & Elmegreenl il998: 
Elmegreen & Elmegreen 2001). In such a hierarchy, the smallest 



observed structures are the youngest ones. This trend is visible 
in Fig. [8] were the spatial distribution of the youngest cluster 
sample shows most clustering and the amount of clustering de- 
creases with the age of the cluster sample. We will quantify this 
more precisely in Sect. [53] with the help of two-point autocorre- 
lation functions. 



5. 1 . Radial distributions 

From the 2-dimensional spatial distributions of SFR, gas sur- 
face densities and star clusters (Fig. [8J we derived radial density 
distributions using 13" wide circular annuli and, in case of the 
clusters, masking out the regions not covered by the WFPC2 
fields (i.e. the grey area in Fig.[8](right)). The SFR and gas aper- 
tures cover the entire field of view, so no masking was neces- 
sary for those. In Fig. [9] we compare the radial number density 
distributions of the different cluster samples directly to the SFR 
surface density distribution. Fig. [10] shows the radial distribu- 
tion of the total gas surface density (including helium and met- 
als), and in Fig. Q~T] we show the radial surface brightness dis- 
tributions of the older stellar populations (no maskin g applied), 
derive d from 2MASS J , H and AT-band ima ging dJarrett et al.l 
120031) . In a recent paper, Schuster et al. (2007) also presented ra- 
dial distributions of SFR and gas density in M 5 1 using 20 cm 
radio continuum data for the SFR density and a combination of 
CO and H I data for the gas density. These distributions can be 
used as an independent check to the data tak en from K07 a nd are 
also shown in Figs. [9] & [TO] We note fhat lSchuster et all d2007l) 
have corrected their distributions (i.e. the galactocentric radii) 
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Fig. 8. Left: Locations of the iKennicutt et all (|2007) 13" apertures used to trace the SFR and gas surface density in M 51. The 
greyscale indicates the mean SFR surface density. The spiral arms from Fig [7] are overplotted for reference. The two dashed circles 
are centred on the nucleus of M 51 and have a radius of 2.5 and 5 kpc, respectively. Right: Distribution of the cluster samples with 
different ages across the disc of M 51, in the same scale as the figure on the left. Filled blue (greyscale: black) circles: clusters with 
log(age/yr)< 7.0. Open green (light grey) circles: clusters with 7.0 < log(age/yr) < 7.5. Purple (dark grey) crosses: clusters with 
7.5 < log(age/yr) < 8.6. The regions not covered by the WFPC2 fields, i.e. for which we have no age estimation of the clusters, are 
indicated in grey. The spiral arms from Fig [7] are overplotted for reference. 



for the inclination angle of the disc of M 51, while the cluster 
densities and the densities from K07 have only been divided by 
an additional factor of 1.07 to correct for the effect of the in- 
clination on the projected area of the apertu re s/annuli. With an 
inclination angle of only 20° dTullv|[T974b . however, this dif- 
ference i^jie^hgjbk for_our purposes. Figures [9] & \10\ show 
that the Schust er et al.l d2007) distributions are smoother than the 
K07 distributions. This is consistent with the radio data having a 
lower resolution (namely a half power beamwidth of 15" for the 
SFR data and 13" for the gas data, compared to the 13" apertures 
of K07). The offset in the gas s urface densities (Fig.fTOli between 
K07 and Schu ster et al.l d2007l) is probably due to differences in 
the conversion factors used to calculate the PL column density 
based on the CO data. 

From Fig. [9] & [10] we see that both the SFR and gas density 
are peaking in the centre, and show two other distinct peaks at 
~2.5 and ~5 kpc. These peaks are not simply a side effect of 
tracing a larger part of the spiral arms at these distances, as they 
are not clearly visible in the /^-band surface brightness distribu- 
tion (Fig. Q~T]i, while this band was used to define the spiral arms 
(Sect. 14. U . The locations of the peaks however, coincide with 
the locations of the four most active star forming regions in the 
spiral arms of M 51, as can be seen in Fig. [8] (left). 

We see from Fig. [9] that the radial distribution of clusters 
younger than 10 Myr (sample 1) follows the SFR density distri- 
bution to a high degree (including the 2 peaks). Sample 2 shows 
a much flatter distribution but is still peaking towards the cen- 
tre, followed by a dip within the inner ~1 kpc, while the old- 
est cluster sample (sample 3) shows a flat distribution. These 
conclusions hold if we would use the Padova isochrones, and it 
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Fig. 9. Radial number density distributions of the different clus- 
ter samples (right axis), compare d to the radial distrib ution 
of the SFR surface density from IKennicutt et alj ((2007) and 
Sch uster et all d2007l) (left axis). 



agrees very well with the results of Hwang & Led (|2008), who 
find peaks at ~3 and ~6 kpc in the radial distribution of their 
"class 2" clusters (i.e. clusters which are more distorted and/or 
have multiple neighbours) and an increasing flattening of the ra- 
dial distribution with increasing cluster colour (i.e. with increas- 
ing age, although some of their reddest clusters might be very 
young and extincted). 
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Fig. 10. Radial distributions of the total gas (i.e. including he- 
lium and m etals) surface density. T he solid curve is derived from 
the data of Kennicutt et al. (2007), while the dashed curve is 
from Fig. 4 of lSchuster et al.l (120071) . 




Fig. 11. Surface brightness distribution of M 51 in the 2MASS 
J, H and K-b&nd. 



The peaks in Esfr, 2 gas and Ectuster at Rc ~ 5 kpc are 
consistent with an enhanced star formation activity caused by 
corotation (Elmegreen et al. 1989), since the corotation radius 
of M 51 is at j? n = 5.4 kpc, assuming a circular velocity of 
V - 200km s _1 dGarcia-Burillo et al.ll 19931) and a pattern speed 
of £2 p = 37km s^'kpc" 1 jzimmer et al.l 20041). The peaks at 
Ra ~ 2.5 kpc could be explained by the presence of the 4:1 
resonance, which is were Q - k/4 equals the pattern speed, CI be- 
ing the angular velocity o f the stars and interstellar material an d 
k the epicyclic frequen cy dContopoulos & G rosb0l 1986, 1988). 
Elmegree n et al.l (1 19891) have found evidence for a spiral arm am- 
plitude minimum at this resonance, consistent with the gaps, or- 
bital loops and stellar orbit crowding, tha t are expected for this 
resonance (Contop oulos & Grosb0llll986l) . It is also suggested 
that orbit crowding and loops result in more gas shocking. This 
would lead to enhanced star formation at the 4: 1 resonance ra- 
dius, resulting in an enhancement of Zsfr at ~2.5 kpc, which 
is then followed by an enhanced E c iuster for the youngest cluster 
sample. 
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Fig. 12. Normalized distribution of the azimuthal distance of the 
clusters to their closest spiral arm. 
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Fig. 13. Normalized distribution of the kinematic ages of the 
clusters from their closest spiral arm, for the clusters with a 
galactocentric distance < 4 kpc. 



5.2. Azimuthal distributions 

If the formation of star clusters is triggered by ga s shocking due 
to the passing spiral density wave (Roberts 1969), we expect the 
youngest star clusters to be located closer to the spiral arms com- 
pared to older clusters. Using the spiral arms from Fig. [7] we de- 
rived for every cluster their azimuthal distance from their closest 
spiral arm, assuming circular orbits. No corrections for the po- 
sition angle (PA) and inclination angle (z) were applied, since 
these will practically be negligible for t he almost fa ce-on orien- 
tation of M 51 (PA = 170°, i = 20°. If^[l974l) . In Fig.QJ 
we show the distribution of the azimuthal distances for the three 
cluster sa mples (sometimes re f erred t o as the "circular distribu- 
tion", e.g. Boeshaar & Hodgi d!977l) ). normalized by the total 
number of clusters in every sample. The figure shows that all 
cluster samples are peaking at the location of the spiral arms 
(i.e. an azimuthal distance of 0°). Although the youngest clusters 
(sample 1) show the most smooth distribution (i.e. with less scat- 
ter) compared to the other two samples, K-S tests do not indicate 
statistically significant differences between the samples. This re- 
sult is robust against the variations of the SSP model parameters, 
used to derive the cluster ages, as described in Sect. 13.21 
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The fact that the oldest cluster sample is peaking at the centre 
of the spiral arms is likely a side effect of the enhanced cluster 
formation around the co-rotation radius (Sect. 15. U . Around the 
co-rotation radius the drift velocity away from the spiral arms 
is very low, meaning that older clusters once formed in these 
regions are still located close to the arms. The azimuthal dis- 
tance of the clusters from the spiral arms expressed in terms of 
"kinematic age" will therefore be a better description of the dif- 
ferences in the circular distributions between the three cluster 
samples. In Fig.[13]we show the circular distributions in units of 
kinematic age, 7\i n , which is the azimuthal distance to the clos- 
est spiral arm, divided by the local drift velocity at the location 
(Rg) of the cluster. A negative T^ R implies that the cluster is 
still approaching its closest spiral arm. In Fig.[T3]we only plot- 
ted clusters with Rq < 4 kpc, i.e. well within the co-rotation 
radius, since close to co-rotation 7\i n will diverge and outside 
the co-rotation radius the "kinks" or near-circularity of the spiral 
arms could make the intersection of the cluster orbit and spiral 
arm ambiguous. Note that this prescription for T^ n does not ac- 
count for the passing of multiple spiral arm, e.g. an old cluster 
which has passed numerous spiral arms in its orbit can still be 
appointed a very low kinematic age, if it currently happens to be 
located close to an arm. 

This notwithstanding, we see in Fig.[l3]that the oldest cluster 
sample is peaking at older kinematic ages (10-15 Myr) com- 
pared to the younger two cluster samples, which peak at 0- 
5 Myr. We also see that sample 1 has a larger fraction of clus- 
ters at the youngest kinematic ages left of the peak compared to 
sample 2. The kinematic ages of the peaks in the distributions 
are younger than the typical ages of the clusters in the samples, 
indicating that the majority of the clusters formed ~ 5-20 Myr 
before their parental gas cloud would have reached the centre 
of the spiral arm. This is in agreement with the sharp edge of 
the dust clouds, seen on the HST/ACS image (see e.g. the dust 
lane covered by field 1 and 2 in Fig. [TJ- For a collapse time 
of the order of 10 years we expect most contracting clouds to 
have forme d a cluster well before the centre of the spiral arm 
is reached. Tambur ro et al. I (120081) find a timescale of 3.3 Myr 
in M 5 1 for the transition from H I emission, via CO emission to 
24/iin emission from current star formation. Since H I first has to 
condense to H2, the collapse time has to be shorter than 3.3 Myr. 

The circular distributions are reasonably symmetrical around 
their peaks, which is consistent with the circular distributions 
found for H II regions in other spiral galaxies (e.g. in NGC 4321, 
lAnderson et al.|[l983l) and it is in qualitative agreement with the 
density wave models which predict a more s moothed su rface 
brightness distribution across spiral arms (e.g. lBashl[l979h . as 
opposed to the models which predict an asymmetric distribution 
dRobertslll969l) . 

5.3. Clustering of star clusters 

The degree of clustering for different cluster samples can be 
quantified and compare d by using the two-point correlation 
function (Peebles 1980, § 45). Applied to clusters from the 
same cluster sample, the two-point correlation function becomes 
an autocorrelation function. I n the fol l owing , we broadly fol- 
low th e method described in Peebles ( 1983) and IZhang et alj 
(2001). Zhang et al. applied two-point correlation functions to 
star clusters in NGC 4038/39 (the "Antennae"). The autocorre- 
lation function is defined as: 
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Fig. 14. The autocorrelation function as function of radius for 
the three cluster samples. 



where n,(r) is the number density of clusters found in an aper- 
ture of radius r centred on, but excluding cluster i. N is the to- 
tal number of clusters and h is the average number density of 
clusters (excluding the grey regions in Fig. [8] (right)). In gen- 
eral, £(r) is defined such that «[1 + ^(r)]dV is the probability 
of finding a neighbouring cluster in a volume d V within a dis - 
tance of r from a random cluster in the sample (Peebles 1980). 
For our spatial distribution of clusters this means that 1 + £(r) is 
a measure for the mean surface density within radius r from a 
cluster, divided by the mean surface density of the total sample 
(i.e. the surface density enhancement within radius r w.r.t. the 
global mean). Therefore, a random distribution of clusters will 
have a constant 1 + £(r) = 1, while for a clustered distribution 
l+£(r)>l. 

Since the star clusters in M 5 1 reside primarily in the disc, 
while the observations cover a larger projected area, all cluster 
samples will show some degree of autocorrelation and the abso- 
lute scaling of 1 + £(r) will be arbitrary. We can, however, com- 
pare the relative scaling of 1 + £(r) between the different cluster 
samples as well as the different functional forms. 

In Fig. [14] we show the autocorrelation functions as a func- 
tion of radius for the three cluster samples. Before the auto- 
correlation was calculated, the original ACS coordinates were 
mapped onto a new image with 1 " resolution to reduce comput- 
ing time. We took into account that as a result of this binning, 
one p ixel could con tain multiple clusters. For the errorbars we 
used (iPeeblesll 19801 § 49): 



6{f) = 



' j N > 

2Z"P (r) 



-1/2 



(2) 



(1) 



where n p (r) is the number of pairs formed with the central clus- 
ter i of the current aperture, and the factor 1 /2 accounts for not 
counting every pair twice. 

Fig. [T4l shows that for r > 4" (~ 160 pc) the clusters with 
ages < 10 Myr (sample 1) are most clustered (largest £(r)) and 
that the amount of clustering decreases with age. This is consis- 
tent with the visual interpretation of Fig. [8] (right). This result is 
robust against variations in the isochrones (Padova or Geneva) or 
metallicity range used in the SSP model fitting of our clusters. 
Below r ~ 10" (i.e. ~ 400 pc) we see a steady increase in the 
clustering of sample 3, surpassing sample 1 below r ~ 4". We 
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also see that below this radius the amount of clustering for sam- 
ple 2 decreases. However, if we would use the Padova isochrones 
we find that the autocorrelation function of sample 3 is much 
flatter, while the turnover of sample 2 decreases. Due to this de- 
generacy, we can not draw any firm conclusions on the age de- 
pendency of the autocorrelation function on scales below ~ 4". 

Moreover, on these smallest scales we can be biased by 
selection effects in the original sa mple selection of resolved 
sources of Scheepmake r et al] d2007l) . In this sample selection, 
sources with a neighbouring source within 5 ACS pixels (i.e. 
0.25") were rejected, possibly lowering the amount of autocor- 
relation in the first radius bin in Fig. [T4j Crowding effects and 
the high background in the most compact regions in which the 
star clusters form, will also introduce a bias against the detection 
of clusters on the smallest scales. This bias will be stronger for 
younger clusters, since in a hierarchical (i.e. "scale free") model 
of star formation the youngest structures are also the most com- 
pact ones. Correcting for such biases will increase the autocor- 
relation function of our youngest sample (sample 1), as we will 
show below. 

The general trend of 1 + £ decreasing with radius is ex- 
pected for hierarchical (i.e. "scale free" or "fractal") star fo rma- 
tion (e.g. iGomez e"taT]|1993t lLarsonll 1995b iBate et al.lll998h . If 
we consider a compact, isolated grouping of newly formed star 
clusters, we expect such grouping to disperse with a typical ve- 
locity dispersion < 10 kms" 1 dWhitmore et alj |2005). For a 10 
and 30 Myr old population of clusters, autocorrelation is there- 
fore expected to increase on scales smaller than roughly 100 and 
300 pc, respectively (which are lower limits, considering the ini- 
tial size > pc of the GMC that formed the grouping). On larger 
scales, the autocorrelation for an isolated grouping of clusters is 
expected to level off. However, this will only hold if the group- 
ing was not part of a larger hierarchy of star formation. In a 
hierarchical model, on the other hand, the clusters will be corre- 
lated with other clusters from other groupings on larger scales, 
together forming a group higher up in the hierarchy of star for- 
mation. The smooth decline of 1 + £ with radius, observed for 
our cluster samples on scales r > 3", is therefore consistent with 
these populations being part of this hierarchy. The velocity dis- 
persion will cause an overall decre ase of the amoun t of cluster- 
ing of the star clusters with time (Ba te et aT1ll998l) . consistent 
with Fig. [H 

For a hierarchical, or fractal, distribution of stars or clus- 
ters, the autocorrelation function will show a pow erlaw depen- 
dency with radius of the form 1 + £(r) oc r n dGomez et alj 
1993). For such a distribution, the total number of clusters (AO 
within an aperture increase with radius as oc ft ■ r 2 = r' ,+2 , 
which shows that r\ is related to the (two-d imensional or pro- 
jected) fractal dimension D2 as D2 — rj + 2 (Mandelbrot 1983; 
iLarsonl [l995h . The distribution of star formation and inter- 
stellar gas over a large range of environments is observed to 
show a (three-dimensional) fractal dimension of D3 ~ 2 . 3 (e.g . 
Elmegreen & Falga rond 1 19961; lElmegreen"& Elmegreen 2001) 
Projecting this three-dimensional fractal on a plane will result 
in a two-dimensional fractal dimension D2 ~ D3 — 1 — 1. 3 
dMandelbrotll983blBeecrJl992blElmegreen & Elmegreenl200ll) . 
In the spiral galaxy NGC 628. lElmegreen etaLl (2006) observe 
star-forming regions with D2 = 1.5. 

If we fit powerlaws to the straight parts of the autocorrelation 
functions in Fig.[l4](4" < r < 50"), we find that the decreasing 
trend of sample 1 and 2 can be approximated with a power law 
index of around -0.4 (fractal D2 = 1.6). Sample 3 has a lower 
index of around -0.8 (D2 = 1.2), but as mentioned before, this 
index depends strongly on the stellar isochrones used: using the 



Padova isochrones, some clusters from sample 2 are fitted with 
older ages (see Fig.|4]i and move to sample 3, flattening the auto- 
correlation function to a power law with an index of around —0.4. 
Intere stingly, an index of -0.4 is similar to the index Zhan g~et~aTI 
d200ll) find for candidate young stars in the Antennae galaxies 
(namely an index of -0.41), while the same authors find steeper 
slopes for the autocorrelation functions of young (5-160 Myr) 
star clusters (namely indices between -0.74 and -1.06). In an 
attempt to minimize the possible bias due to crowding of the 
youngest clusters, we also determined the autocorrelation func- 
tion of all clusters in sample 1 with log(M/M ) > 3.7. This 
sample will have a higher completeness in the densest star form- 
ing regions. For this mass-limited sample we find a power law 
index of -0.8 (D 2 = 1. 2, Geneva) or -0.7 (D 2 = 1.3, Padova), 
similar to the indices o f lZhangetalJd2001h and consistent with 
the fractal dimension of ~ 1.3, mentioned above. This suggest 
that the flatter slope of the full sample 1 population of clusters 
is indeed biased, and that the unbiased index of the autocorrela- 
tion function of young clusters is similar for different galaxies. 
In other words, this hints towards a universal fractal dimension 
of hierarchical star formation. More studies in different galaxies 
and of different nature are necessary to strengthen this point. 



5.4. Spatial correlation between star clusters and star 
formation 

The crosscorrelation function (Eq. [TJ can also be used to study 
the correlation between star clusters and the flux from any image. 
We can rewrite Eq.[TJas: 



1+ ^ = ^5>)' 



(3) 



following Zha ng et all d2001l their Eq. 2), where /j(r) is now the 
intensity (i.e. flux per pixel) in an aperture with radius r cen- 
tred on cluster i, and / is the mean intensity over the whole im- 
age. For both the calculation of /;(/*) and / the areas not cov- 
ered by {/-band imaging are masked out. To study the cross- 
correlation between star clusters and st ar formation, we use d 
the 20 cm radio continuum (RC) map dPatrikeev et al. 2006), 
and t he HST/ACS Ha mosaic image of M51 dMutchler et al.1 
120051) . which we continuum subtracted using a combination of 
the ACS Vf s55w and If%uw imag es. The 20 cm RC map was 
also used bv lSchuster et al.l ((2007) to derive the radial SFR dis- 
tribution (shown in F ig. [9]). The SFRs from the apertures of 
iKennicutt et al.l d2007l) (Fig. [8} are not suitable to be directly 
used in Eq.|3]due to their incomplete spatial sampling. 

The 20 cm RC is tracing star formation through the syn- 
chrotron radiation coming from cosmic rays, emitted by super- 
novae explosions of massive stars. The exact relation between 
SFR and RC emission, however, is not well defined. It depends 
on the scaling between RC and far-infrared dust continuum 
emission, of which the gene ral form and radial dependencies 
are st i ll under debate (see e.g . lHelou et al.ll9 85; Mars h & Heloul 
1991 iNiklas & Beckl fT^97t iNiklasI 1 1 997t iMurgia et all 12005b 



Mur phy et al.ll2.Q06l 120081) . For our current study, however, we 



are not interested in the exact scaling between RC or Ha flux 
and SFR. We simply use the RC and Ha flux to study how star 
clusters and two spatially complete tracers of star formation cor- 
relate for differently aged cluster samples. 

Before we applied the crosscorrelation function (Eq. [3]) to 
the 20 cm RC image (shown in Fig. [TBI), the image was cropped 
to match the ACS mosaic frame and resampled to a pixel scale 
of 1" (the original image had a pixel scale of 2" and a beam 
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Fig. 15. Radio continuum (20 cm) image of M 51, with the spiral 
arms of Fig. |7]overplotted for reference. 



size of 15")- Due to computational issues the Ha image, with 
a resolution of 0.05" was binned by a factor of 10 to a pixel 
scale of 0-5". For t he uncertainty in 1 + £(r) we again followed 
Zhan g et all ((2001), adopting only the statistical uncertainty in 
the number of clusters per sample (AO, leading to: 



S(r) = AT 1/2 (1 +£(r)). 



(4) 



In Fig. [16] we show the resulting crosscorrelation functions 
between the three cluster samples and the RC flux. The crosscor- 
relation functions are flat on scales < 10", which is very likely 
a result of the beam size of the radio data. It is evident that the 
youngest cluster sample is most correlated (i.e. largest ^(r)) with 
the RC emission on all scales smaller than ~ 100" (i.e. ~ 4 kpc), 
and that the correlation drops with the age of the cluster popu- 
lation. A random distribution of clusters would have a constant 
£(r)+l = 1 . In other words, most of the radio flux is concentrated 
around the youngest clusters, as expected. 

The crosscorrelation function between the 20 cm RC emis- 
sion and the youngest cluster sample drops to half its peak value 
at R ss 2 kpc. This shows that the RC emission is tracing star for- 
mation in a very diffuse way, consistent with the long path length 
of the cosmic rays, before they stop emitting RC emission. 

In Fig. [17] we show the crosscorrelation functions between 
the three cluster samples and the (continuum subtracted) Ha 
flux. The youngest clusters are strongly correlated to the Ha flux. 
The correlation peaks at the resolution limit of R = 0.5" ~ 20 pc 
and drops to half the peak value at R » 100 pc. This radius is 
consistent with the typical expans ion speed (~10 kins" 1 ) of the 
ionized bubbles of hot stars (e.g. La mers & CassineIlHll999l p. 
369). At a radius larger than the typical travel distance of an Ha 
bubble in 10 Myr (i.e. ~100 pc) the correlation is expected to 
level off. For sample 2 clusters the correlation with Ha is con- 
siderably less and the sample 3 clusters shows a crosscorrela- 
tion function with Ha consistent with a random distribution (i.e. 
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Fig. 16. Crosscorrelation function between the three cluster sam- 
ples and the 20 cm radio continuum image. 
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Fig. 17. Crosscorrelation function between the three cluster sam- 
ples and the (continuum subtracted) Ha image. Note the dif- 
ference in X- and Y-scaling compared to Figs. [16]&[T8] On a 
larger scale, the three curves would converge to £ + 1 = 1.1 at 
R « 120" « 5 kpc, similarly to Fig. [TBI 



£(r) + 1 = 1, although a slight offset from 1 is expected due the 
large scale structure or M 51, see below). 

As a consistency check, we also derived the crosscorrela- 
tion functions between the star clusters and the flux from the 
2MASS //-band image, which we used to define the spiral arms 
(Fig. [7]l. Since this near-infrared data is mostly tracing older stel- 
lar populations, no enhanced correlation with young star clusters 
is expected. In Fig.[18]the crosscorrelation functions are indeed 
shown to be indistinguishable. All three clusters samples, how- 
ever, show an increased correlation towards smaller radii. This 
is likely an effect of the large scale structure of M 5 1 . Both the 
cluster and //-band flux density are higher in the disc and spiral 
arms of M 5 1 , while the observations (and thus the normalisation 
factor / in Eq. [3]l cover a larger area. This leads to some degree 
of correlation towards smaller scales. The fact that this correla- 
tion is similar for all cluster ages confirms that the //-band flux is 
a good tracer of stellar density and not of current star formation. 
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Fig. 18. Crosscorrelation function between the three cluster sam- 
ples and the 2MASS //-band image of M 5 1 . 



6. Summary and conclusions 

We have studied the spatial correlations between star clusters 
and star formation across the spiral galaxy M 5 1 . We combined 
published data of the resol ved star formation law in M 5 1 (SFR 
and gas surfa ce densities, [Ke nnicutt et al. 2007) with the star 
cluster data of IScheepmaker et al.l (120071) . We combined the star 
cluster data with new HST/WFPC2 U F3i6w data of 1580 clus- 
ters in order to estimate ages, masses and extinctions of these 
clusters using SSP models. We used radial and azimuthal distri- 
butions and cross-correlation functions to study the spatial cor- 
relation between star formation, star clusters and spiral structure. 
We used auto-correlation functions to study how the clustering 
of star clusters depends on age and spatial scale. The main con- 
clusions of this work can be summarized as follows: 



1. We find that age determinations of young (< 100 Myr) 
star clusters based on broad-band UBVI photometry can de- 
pend strongly on the adopted stellar isochrones used in the 
SSP models. Using models based on the Padova isochrones, 
we generally find older, and therefore more massive, clus- 
ters compared to the same models based on the Geneva 
isochrones (Fig.|4]i. 

2. The radially averaged distributions of star formation and gas 
surface density in M 5 1 peak in the centre and subsequently 
at ~ 2.5 and ~ 5 kpc. These distributions are matched by 
the number surface density distribution of our youngest star 
cluster sample (ages < 10 Myr, Fig. [5), while for the older 
clusters the distributions are more flat. We note that the loca- 
tion of the peaks coincide approximately with the locations 
where the spiral arms meet the 4:1 resonance and the co- 
rotation radius. 

3. The distributions of azimuthal distances from the spiral arms 
show that most of the young and old star clusters are lo- 
cated at the centre of the spiral arms, with the youngest star 
clusters showing the most smooth distribution (Fig. [121 ). In 
terms of kinematic age and for Rq < 4 kpc, the azimuthal 
distribution of the oldest clusters peaks at the largest separa- 
tion downstream of the spiral arms (10-15 Myr), while the 
majority of the younger clusters are within 0-5 Myr down- 
stream of the arms. This indicates that the majority of the 
clusters formed ~ 5-20 Myr before their parental gas cloud 
would have reached the centre of the spiral arm. 



4. We find that the amount of clustering of star clusters de- 
creases both with spatial scale and age (Fig. RBI, consistent 
with hierarchical star formation. The slope of the autocor- 
relation functions of the youngest clusters (log(age) < 7.5) 
is consistent with a projected fractal dimension of 1.6. We 
find a lower fractal dimension of 1.2 for the clusters with 
log(age) < 7.0, if we use a higher mass limit and thus, for 
this smaller mass range, a more complete sample. On spatial 
scales smaller than 400 pc, the amount of clustering of the 
oldest cluster samples is very uncertain due to uncertainties 
in the assumed SSP model parameters. 

5. The locations of the youngest star clusters are highly cor- 
related with the flux from 20 cm radio-continuum emission 
and continuum subtracted Ha emission, but not with near- 
infrared //-band emission. The correlation with the radio 
emission is more diffuse compared to the correlation with 
Ha, and both correlations decreases with the age of the clus- 
ter population. 



In a follow-up paper (Scheepmaker et al. 2009) we will focus 
on the quantitative comparison between star and cluster forma- 
tion, discussing the fraction of total star formation taking place 
in optically visible star clusters. 
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